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ABSTRACT 

Growing interest in the recovery of astrovehicles and artificial 
planets has led to investigations of surface cooling at hypervelocities 
(Mach numbers greater than 15). The simulation of these problems 
using models in electromagnetic shock tubes has offered a direct 
approach toward verification of predicted heat transfer rates. This 
thesis is an analysis of a proposed experimental apparatus which will 
measure the stagnation point heat transfer rate to models in highly 
ionized flows in the electromagnetic shock tube. 

It is proposed that stagnation point temperature measurements be 
made with a photoconductive cell by measuring the infrared radiation 
emitted from a metallic (gold) film placed on a quartz window in the 
nose of the model. From the measured variation of temperature with 
time, the heat transfer rate can be determined. The transmission 
path of the infrared radiation is along a sapphire light pipe to a lead 
sulphide detector external to the model. The device is capable of 
operating under conditions where high ionization levels make the use 
of resistance film gages impractical, and it may also be operated in 
the presence of varying magnetic fields. The device is capable of 
detecting heat transfer rates between approximately 2000 and 
55,000 ane, The lower value is determined by the available 
period of observation (40 microseconds) while the upper value is 
limited by the detector cell time response, melting temperature of 
the gold film and the available period of observation. The response 
time of the system limits observation periods to those equal to or 


greater than the available minimum detector cell time constant (about 


ii 





10 microseconds). Expected heat transfer rates are calculated and an 
optimum operating wavelength is derived. Component specifications 
are determined and calibration techniques are suggested. 

The system concept and preliminary design work were originated 
by Dr. R. W. Ziemer at the Physical Research Laboratory, Space 
Technology Laboratories, Inc., Los Angeles, California. The system 
is presently being constructed and is expected to be in operation some- 
time in 1959. This analysis was made by Larry G. Valade, Lieutenant, 
U.S. Navy, and William Rees Phillips, Lieutenant, U.S. Navy, at the 
Physical Research Laboratory through the kindness of Dr. Milton U. 


Clauser, Director. 
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1. Introduction. 

With the current interest in space vehicles and re-entry heating, 
the problem of measuring the heat transfer rate to a blunt bodyina 
highly ionized, magnetoaerodynamic flow has become increasingly 
important. A device is proposed for the measurement of the stagna- 
tion point heat transfer rate on a test model and was designed for use 
in the electromagnetic shock tube shown in Fig. 1. The basic concept 
consists of placing a thin metallic film on the nose of the model, shown 
in the shock tube in Fig. 2, and then observing the infrared radiation 
at a selected wave length with a photoconductive detector cell. The 
ultimate objective is to verify the magnetic cooling effect predicted 
by Meyer [1] and Bush [2]. This design is a continuation of the mag~ 
netoaerodynamic research work performed by Dr. R. W. Ziemer and 
W.B. Bush, Space Technology Laboratories, Los Angeles, California, 
in determining the magnetic field effects on hypervelocity flow about a 
blunt body (3). 

Resistance film heat transfer gages have been utilized with a 
great deal of success for this type of measurement [4]. However, at 
the higher Mach numbers (above approximately 15), the ionized parti- 
cles striking the resistance gage generate noise greater than the 
normal signal, making the system impractical for such flows. 

Successful observations of explosion temperatures by photo- 
multiplier color pyrometers (5) Suggested radiometric measurement 
of the surface temperature of projectiles fired at hypervelocities. 
This method is being employed at the U.S. Naval Ordnance Test 
Station. Cline taake — Galitornia: [6] utilizing electron photomultiplier 


tubes operating at wave lengths of 0.55 to 0.75 micron. The application 





of this technique in the infrared spectrum to the magnetoaerodynamic 
heat transfer problem was suggested by Ziemer. 

The problem of detecting infrared radiation from the modelina 
shock tube with an approximate stream temperature of 8000°K, and 
shock layer temperatures up to 35,000°K, presents many interesting 
considerations. In addition, the short duration of steady state condi- 
tions (15 to 40 microseconds), the necessity of shielding the detector 
cell from the free stream and shock layer radiation, and the need for 
electromagnetic shielding of the electronic components, places stringent 
requirements on the system. 

Initial calculations indicated that model nose film temperatures 
from 300 to 1200°K could be expected. The wave length vol pear 
radiationat these temperatures lies betweentwo andten microns. The 
greater sensitivity of lead sulphide cells (over other photoconductive 
detectors) will be required for maximum system response. Typical 
spectral response characteristics|7| of a lead sulphide cell are shown 
in Fig. 3. Peak response of this cell lies between 1.9 and 3.5 microns. 
The high temperature radiation from the free stream gases(8| exhibits 
an apparent minimum at a wave lengthof 3.1 microns. These considera- 
tions suggest an initial spectrum of interest between wave lengths of 
two and four microns. 

In designing the system, the following will be considered, 

(1) A study of the shock tube characteristics and evaluation of the 

expected heat transfer rates 

(2) Selection of the metallic film material and thickness 

(3) Determination of a suitable energy transmission path 


(4) Selection of an optimum wave length 





(5) Detector cell performance 
(6) Instrumentation 


(7) Calibration. 





2. Conditions of Heat Transfer Measurements in the Electromagnetic 

Shock Tube. 

In order to design a heat transfer measurement system it will be 
necessary to first examine those conditions under which the measure- 
ments will be made. It will be shown later that the system must be 
capable of measuring extremely high transfer rates up to 35,000 Btu/ 
2 aes The measurements must be made within the 15 to 40 micro- 
second steady-state times inherent in electromagnetic shock tubes. The 
device must also be insensitive to strong electromagnetic radiation 
from the shock tube gases and the electrical circuitry. 

The test model is shown in Fig. 2. The physical dimensions 
closely approximate those of models used in previous work with the 
shock tube. The operation of this particular shock tube has been 
described in detail by Ziemer. [9] The following basic principles are 
outlined for the reader unfamiliar with this device. 

A high voltage electrical discharge at the base of the tube suddenly 
adds driving energy to a portion of the gas. As the gas rapidly expands, 
a simultaneous magnetic pinch forces the gas down the shock tube. A 
very strong shock wave with Mach numbers of 40 or more is generated. 
The main equipment associated with the tube consists of the shock tube 
assembly, electrical spark gap switch, high voltage power supply, 
condenser bank, trigger pulse circuitry, high current pulser for 
applying the magnetic field at the nose of the model, and a vacuum 
system for evacuating the tube. The complete system is shown 
schematically in Fig. 4. 

Various shock tube operating conditions are obtained by evacuating 


the tube to various initial pressures in the range from 30 to 300microns. 
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the tube to various initial pressures in the range from 30 to 300 microns. 





The tube is usually operated with air, but could be filled with any 
gas of interest. Shock Mach numbers are approximately inversely 
proportional to the initial tube pressure at the time of firing. 

In order to fire the tube, a high voltage supply charges the capaci- 
tor bank to about 25 kilovolts, which is applied between the ring and 
lower spark gap electrodes. The spark gap prevents current flow until 
the 15-kilovolt pulse is applied to the lower switch plate. Then the 
voltage between the lower switch plate and the trigger wire exceeds the 
breakdown potential, and this initial discharge ionizes the gas in the 
spark gap, which becomes conductive. This permits the capacitors 
to discharge, and current flows from the center electrode to the ring 
electrode and down the six copper straps placed around the outside of 
the tube's conical section. The resulting high current produces a 
magnetic field which pinches and further heats the ionized gas in the 
cone, driving the shock wave down the tube. 

A high current can be pulsed through the coil located in the nose of 
the model, resulting in maximum magnetic field strengths of about 
40 kilogauss at the stagnation point. The current is supplied by an 
L-C circuit whose resonant frequency is approximately 2000 cps. 

The shock tube is fired just prior to the time when the field strength 
is at its peak value. This provides an essentially constant magnetic 
field over the 40-microsecond test period due to the relatively long 
500-microsecond period of the L-C pulsing circuit. 

In order to determine the thermodynamic conditions existing in 
the tube, the shock wave velocity is measured experimentally by 


means of the rotating mirror camera. [10] Knowing the shock velocity 





together with the initial pressure or density, is sufficient to calculate 


the characteristic properties of the flow. 


The shock tube flow characteristics, free stream density ratio, 
as functions of shock Mach 


P>/ Po: and the free stream velocity, Us: 
number are determined as follows: 


From the continuity equation, 
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Utilizing the Hugoniot shock wave relationships and the thermodynamic 
properties of high temperature air reported by Gilmore, (1 1} the desired 
quantities are readily obtained for various shock tube initial conditions, 
and these are presented graphically in Fig. 6. 

The operating characteristics of the shock tube using air may be 


summarized as follows: 


Shock Mach number, M. 15 to 40 

Initial density ratio, p/p, eeilne a) 5 10 - 
Free stream density ratio, pa ee 3,60 110) teno eoeiee 
Free stream Mach number, M, 2. 91tOr 3.9 

Free stream velocity, u, 15,000 to 41,000 ft/sec 
Free stream temperature, T, 6000 to 11,000°K 
Shock layer temperature, T; Z0,000 to 35,0007 Kk 
Test duration 15 to 40 microseconds 


Pressure shock tube experiments for determining the stagnation 
point heat transfer rate on similar test models have been reported by 
Detra[(i2] of the Avco Research Laboratory, Everett, Massachusetts. 
It was found that at test shock Mach numbers of up to 17, the following 


empirical relationship described the observed heat transfer rates. 


lay = h. 
q. = = 2000 ,,/ so ae | att -sec 
Ss w300 





where R is the body nose radius, ue is taken as a reference velocity, 
26,000 ft/sec, and hi and he the wall and stagnation enthalpies, 
respectively. 


Making the assumptions that, 


ee oe 
h. > A300 


the expression reduces to 


aN) Pe — — nw ee eaee 
sn 0 R Pee 


Substitution of the previously determined thermodynamic properiic.e 
into the equation above results in the heat transfer rate curve asa 
function of shock Mach number and is shown in Fig. 7. The validity 
of this relationship becomes suspect above the range of Detra's 
experimental verification, because it would seem that the radiation 
effects from the shock layer and free stream gases may be contri- 
buting substantially to the heat transfer rate. These effects arc 
probably just becoming significant at the highest velocities observed 
by Detra. 

It should be noted that the Detra equation is valid for atmospheric 
free stream or actual flight conditions and substituting values obtained 
in the electromagnetic shock tube into this expression will yield errone- 
ous values. Free stream temperatures in the shock tube are about 
8000°K, while actual flight stream temperatures will be of the order of 


100°K and will vary with altitude density. When this was taken into 





account, preliminary calculations indicated that the heat transfer rate 
in the shock tube will be approximately 25 per cent greater than that 
shown in Fig. 7. 

An analytical study by Meyerott (13] reported the effect of radiation 


heat transfer is governed by the relation 


qu 3/2 eae 


while the convective heat transfer appears as 


qo oif2 n Ve we 


The expression for convective heat transfer agrees closely with 
Detra's expression. However, a comparison of the two heat transfer 
rates, qd. and qa? within the range of interest indicates that the 
radiative contribution is of the orderof ten per cent or higher of the 
convective portion. This means that the estimated heat transfer rate 
in the shock tube will be greater than that calculated for pure con- 
vective heating. Hence, due to both the "artificial" conditions in the 
shock tube, and the radiative heat effects, the total heat transfer rate 
to the test body in the shock tube will be appreciably greater than that 


predicted by Detra's equation. 





57 ) Lbing PF lime Requirements. 

Several factors must be considered in selecting the metallic film 
to be placed on the test body. A maximum temperature rise is desired 
in the short time available for measurement in order to obtain a high 
emissive power. This requires a film material of high thermal diffusi- 
vity, evaporated on a substrate of low thermal diffusivity. Moreover, 
the film should be of minimum thickness so that the temperature 
gradient within the film will be small. 

The emitted radiation from the film must be observed from behind 
the film so that the observed radiation does not have to pass through 
the hot and also radiating gases in the free stream and shock layer. 
However, thin metallic films may permit appreciable transmission of 
radiant energy [14] in the spectrum of interest. In order to have the 
film opaque to this radiation, the film should have a high reflectivity, 

a high absorption coefficient and a certain minimum thickness. There- 
fore, because of conflicting thermal and optical requirements, a 
compromise film thickness must be determined. 

Thermal and optical properties of various film and substrate 
materials must be considered. The desirability of a gold film anda 
quartz substrate is evident from examination of the properties of the 
most promising materials listed below 

The chemical inertness of gold is another factor that makes its use 
attractive. Certain grades of quartz are opaque above wave lengths of 
2.5 microns, but selected fuzed quartz has excellent transmission out 
to about 4 microns. (15) 

These properties and those physical constants used throughout the 


thesis were largely obtained from American Institute of Physics 
Handbook, McGraw-Hill, New York, 1957. 
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Thermal Specie Spectral Absorption 


Material Diffusivity Heat Emissivity Coefficient 
£t“/sec Btu/hr -ft-°F 274 2-4 Sua 
Sra ee se ey, © 0.031 0.03 20.0 
Siete © eaieeconn 0.032 0.03 3.4 
Quartz Ns Wee 0.18 -- -- 
Pyrex 7.00 x 10 ° 0.20 -- Ss 
Sapphire a aM se UO = Oate 0.03 -- 


The transmissivity of the gold film must be investigated in order 
to determine an optimum thickness which will permit a minimum 
temperature gradient from the outer to the inner film face, anda 
maximum ratio of emitted film radiation to transmitted hot gas radia- 
tion. The transmitted energy is made up of radiation from both the 
free stream and shock layer gases and may be expressed in terms of 
the incident radiant energy on the gold film as 


d 
We = (W, + wi)Q shee scp [-21a a 


Where R is the reflectivity constant, a, the absorption coefficient 
and d, the film thickness. 
The ratio of the emitted energy to the transmitted energy for all 


wave lengths may be expressed in terms of Planck's Law as, 
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Assuming black body radiation from the gases, 


Substituting the value of 
2 = ©).03 ie (eye Zi << y ot 


and then evaluating the energy ratio as a function of \ for various 


thickness with 


T, = 35,000°K 

T, = oOOUOt kk 

Ree Ogos 

a = f(\) from Fig. 8 
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yields the curves of Fig. 9. This graph indicates that a thickness of 
5000 A gives a ten to one ratio of emitted energy to transmitted energy 
for wave lengths of interest. 

The temperature profile in the substrate for a 5000 A thick film 
is determined next as a function of time and heat transfer rate. The 
differential equations describing the heat flow are well known and the 
solution has appeared in several references. (16] : {1 7] The method of 
solving the problem is to first assume the thin film to be a semi- 
infinite slab as shown in Fig. 10. The solution to this configuration is 
later used in the exact solution for a film of finite thickness ona 
substrate of different eal properties. The differential equation 


describing this temperature rise-time relationship in the semi-infinite 





slab is 
z 
ze = K ouee : for Coc 50} 
ons 
where 
K=— 
p p 
A(t) = Tix t) — Peco) 
and 
T(x, t) = temperature at any point, x, inthe slab at time, t 


The boundary conditions to be satisfied are 


Gike, o) eo (1) 


lime O (x, t) = 0 (2) 
x > 0O 
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Using Laplace transform methods, 


L(8) = @(s) 


then 
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Applying boundary conditions (2) and (3) to evaluate constants Cy and 





C5, it is found that 
Cy = +0 
and 
de --2(2)--,\8 
dx K \s 2\s 
x= 0 
giving 
_ q@iNK  -3/2 
e = Ss 
Zz K 


Taking the inverse transform, [18] the solution is 
NK t x? x 1 xal 
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a en KT AK Z AK t 
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This temperature rise-time versus depth relationship is shown 
graphically in Fig. 11 for the general case, and the surface tempera- 
ture solution is plotted in Fig. 12. If the properties of the substrate 
are used in this solution, the approximate solution is obtained for a 
very thin film on a semi-infinite substrate, for which the presence 
of the film has negligible effect. 

However, the solution to the particular problem of a thin film of 
finite thickness mounted on a semi-infinite slab, as shownin Fig. 13, 
remains to be determined. One solution to this temperature rise-time 


relationship tl 9) 1s 
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where 
W(X) > t) = temperature rise at any point, 
xX): in the film at time t 
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It can be shown algebraically that 
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where 9(x,t) is the previous solution for the semi-infinite slab of 
the film material. 
For evaluation of the interface temperature rise as a function of 
time 
x, = £ 


then, simplifying the general expression yields 


Wl, t) = pe? (1 +b) @[ (2n — Lye, ] 


iM 8 


] 


The outer face temperature rise may be expressed as 
arn 
wy (0,t) = 8(0,t) +24 5 b Ol(ené, t) 
n=l 
Evaluating the general expression for the 5000 A thick gold film 
on a quartz substrate for various heat transfer rates, and assuming 


a room temperature of 300°K, gives the interface temperature versus 


time curves of Fig. 14. 


16 





4. Method of Radiation Transmission. 

The problem of transmitting the radiated energy from the gold 
film to the infrared detector cell will be considered next. Knowing 
the characteristics of a typical detector cell, (7} the minimum detect- 
able power level which must fall on the sensitive area of the detector 
can be determined. The temperature-time history of the film inter- 
face is known for the shock tube operating conditions, and therefore 
the power emitted by the film through the quartz substrate is also 
known. 


The transmission problem, then, is to determine whether 


> 
uC - a Pdetectable 
in which Ws = emissive power from the gold film 
qT; = transmission efficiency of the system 
Zee ae = minimum power detectable by infrared cell 


The problem is reduced to choosing a method of transmission with 


sufficient efficiency such that the lower limit established by 


P 
oy = detectable 


W 
& 


is exceeded or that qT, closely approaches 100 per cent. 

The efficiency of any transmission system will be a function of the 
wave length of the transmitted energy. The operating wave length 
remains to be selected at a later stage of the analysis, and, therefore, 
no valid approximation is available for a probable value of qT; for va 
particular method. Therefore, the optimum transmission efficiency 
of a number of possible transmission methods will be examined, making 


the problem one of choosing the most efficient technique. 
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The first transmission system considered is an optical lens 
arrangement which will focus the emitted energy and direct it to the 
detector's sensitive area. Various lens-mirror arrangements were 
considered and it was found that the maximum transmission efficiency 
of any optical system is about 13 per cent. The geometry of the aero- 
dynamic body restricts the efficiency of this method by limiting the 
collecting area of the first lens in the system. 

The successful use of a sapphire rod as a light pipe, in silicon 
crystal growth control, (20) suggested a similar application to this 
problem. A study of the optical properties of crystals in the infrared 
region indicated that the transmissivity of sapphire, in thin samples, 
approaches 90 per cent for wave lengths up to approximately five 
microns. The proposed transmission system incorporating a sapphire 
light pipe is shown in Fig. 15. Based on this diagram, the trans- 
mission efficiency of the rod may be expressed as 


n/2 


i = 2 cos B exp [-aL Sec i] di 


in which L = length of the rod. From Snell's Law for the refraction 


of light 


sin U2 


Sin 1 My 
in which n is the spectral refractive index of the transmitting 


material. Then 


Z 
| tha exp |-a Lseci] di 
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A lower bound for the transmission efficiency can be determined 


by evaluating the expression 





Z 


ul 
Tm (2 | ano A exp [-a Lsecil di 
l 


where the angle, Los is that angle below which total reflection occurs. 

The spectral coefficient, a, for sapphire appears in Fig. 16. 
The integral above was solved graphically over a range of wave lengths 
of two to five microns, and the results plotted in Fig. 17. 

It is interesting to note that the use of the quartz substrate 
improved the transmission efficiency of the system. The refractive 
index of quartz is sufficiently less than that for sapphire so as to pro- 
vide a focusing effect that permits more energy to fall on the face 
between the substrate and the light tube. 

This study showed that a sapphire light pipe is a highly efficient 
means of transmitting energy in the infrared spectrum. The expected 
transmission efficiency of this configuration can be as high as 95 per 
cent in the near infrared and hence, the light pipe was selected in 


preference to the lens system. 
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5. Comparison of Infrared Signal and Extraneous Radiation. 

The intensity of the signal power from the gold film received at 
the cell must be compared to the incident power emitted by all other 
sources. At the lower wave lengths, considerable power is trans- 
mitted through the film from the shock tube gases. At the longer 
wave lengths, power comparable to the signal power is emitted by the 
light pipe, substrate, and the test model walls. Hence, an optimum 
operating wave length will be determined to maximize the ratio of the 
Signal power to the extraneous power. This investigation will require 
selection of the physical dimensions of the apparatus as it is mounted 
on the shock tube. 

The diameter of the sapphire light pipe should be a maximum, in 
order to achieve greater transmission efficiency. However, this 
dimension is limited by the available space within the magnetic coil 
core. Further a small diameter is desirable to minimize the rod 
surface area, and, therefore, reduce the incident transverse power 
into the rod. Furthermore, it is advantageous to have the gold film 
diameter the same as the light pipe diameter. However, a minimum 
gold film area is required in order to approximate closely the stagna- 
tion point conditions. As acompromise, a diameter of 0.15 inch 
was selected for the initial assembly and subsequent calculations will 
be based on the resulting area for the gold film, substrate, sapphire 
rod, and infrared components. 

The detail design of the test apparatus is shown in Fig. 18. The 
housing above the model will be utilized for the detector cell, circuitry, 
and other necessary electronic components. The housing is made of 


thick copper walls in order to shield the electrical components of the 
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system from stray electromagnetic radiation. The 3/1l6-inch minimum 
thickness used will permit penetration of only two per cent of 5000-cps 
radiation. 

Returning to the problem of maximizing the signal to background 
power ratio, it will be necessary to investigate allpossible sources of 
infrared radiation that will be incident on the cell, and if possible, to 
select a restricted spectral band in which the ratio of the emitted power 
from the gold film to the emitted power of all other sources will be a 
maximum. An interference filter will be required to limit the cell 
spectral response to a discrete bandwidth, permitting an optimum 
signal to background power ratio. Possible sources of undesirable 
noise radiation will include, 

(1) Substrate 

(2) Interference filter 

(3) Incident heat flux on the model's structural surfaces, other 

than at the stagnation point 

(4) Eddy current heating of the film by the magnetic field 

(5) Light tube 

(6) Free stream gases 

(7) Shock layer gases 

sources 1] and 2 are negligible when the magnitude of emissivity 
1s considered for thin layers of quartz and germanium. The radiation 
level from these sources remains at room temperature throughout the 
test period as indicated by Fig. 11. 

The structural source level 3, is small, limited by the shielding 


effect of the copper and will be minimized by highly polishing the tube 
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surrounding the sapphire rod. Also the outer surface of the model 
will be fabricated with an opaque plastic material. 

The heating of the gold film by the eddy currents induced by the 
magnetic field, is investigated and the magnitude of the temperature 
rise determined. The physical picture of the effect is illustrated in 
Fig. 19. The energy dissipated in the film can be determined by 


assuming a uniform flux density over the entire film of 
® = AB. sin wt 


and using the basic relationships 


_ de 
Sana 

2 
a 
ap = = 
Epes 
DS 


From the geometry of Fig. 19, the dissipative power is 


ny Z 
aia. eeae ee dt w" Be 
m m 


4 
_ ee) + 
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evaluated for 
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p = 2.42 pohm - cm 
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yielded 
Z 
P = 0.19 cos wt watts 


The temperature rise in the film was estimated by means of the 


relation 


c,pVAT = P dt 


0 


where t, = 600 uw sec, which resulted in 


AT = 4.1 x Tome °K/cycle 


which can be ignored. 

The sapphire rod, while nearly transparent to infrared radiation, 
has an emissivity comparable to that of the gold film in the spectrum 
of interest. This room temperature radiation from the sapphire rod 
is significant at the longer wave lengths and must be considered in 
choosing an optimum operating wave length. Using the familiar 


expression, 


Tears «rod ui ood ie 
and assuming 
God ao 
then, 
“rod 2 ry 
where 
“rod 7 F(A) 
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from which the emitted power of the sapphire rod at 300°K may be 


expressed as 


+) 
dx(l- Toa) 


W = De eek 
i @ | loi fa 
*P 1 300kx 


rod 


My 


Radiation from the sources 6 and 7 will be transmitted through 
the film as a function of film thickness and wave length as discussed 
in Sect. 3. This energy will be significant at the shorter wave lengths 
and, therefore, should be considered in arriving at hi optimum wave 


length for system operation. For a film thickness of 5000 A, and, 


assuming 


ye ta = €i2 = ] 
qT; = 8000°K 

qT, = 20,000°K 

Gs fa) 

R = 0.983 


the radiation from the gases which is transmitted through the gold film 


can be expressed by 


1.) 
: 2 dy 
Wie =. Sopoieuee mae S ch_]_] 
©XP | 8000kX 
My 
Me | 
a =. ee Cee 0 -b-40) |-27a cy 
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The selection of an optimum wave length will depend on evaluating 


the ratio 


W (A, T) 


W,5(\) + W r(x, T) 


OX) = 


rod 
at different wave lengths and various temperatures, then determining 


that wave length which gives 


r(r\, T) = Ee 


The results of this computation clearly indicate, in Fig. 20, an 


optimum signal at 


KH = 3p 


The selection of a proper bandwidth is made by comparing the 
spectral response of the detector in Fig. 3 with the curves of Fig. 20. 
A wide bandwidth is desirable in order to have maximum incident power 
on the cell. However, a narrow bandwidth is required to obtain a high 
quality signal so that the ratio, r, is greater than some arbitrary 
value corresponding to a satisfactory signal-to-noise ratio. 


It seems appropriate to select a bandwidth 
An = 0.4y 


centered at 
hh = 2.8yp 


in order to take advantage of maximum detector response. An inter- 
ference filter, with characteristics as shown in Fig. 21, will be placed 
in front of the detector cell. This will limit the shorter wave length 


response, andthe detector cell will limitthe longer wave lengthresponse. 
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The power ratio, 


is evaluated at 
Wc ot 
as a function of the gold film temperature and plotted on Fig. 22. This 


indicates that a film temperature of 380°K will be required in order to 


achieve a signal with 
nb 


The incident signal power on the detector from the gold film may 
be calculated from the temperature-time relationships in Fig. 14 by 


the following approximation: 


€ Deaoe th 


Se eee 


We § ee ch 
on h exp RET ] 
detector 


where 


rod | eee 


= Casi 
2 
A -= O21 52cm 
a = 0.8 
Weer 
KH = 2.8yp 


This equation is evaluated for three different heat transfer rates 


and the results appear graphically in Fig. 23. 
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Steady-state flow conditions are reached in the shock tube about 
five microseconds after the tube is fired. In order that the emissive 
power at this time be detected, the cell must be capable of sensing a 
minimum incident signal power of 0.015 microwatts for low shock 


Mach numbers. 
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6. Detector Cell Performance and Instrumentation. 

Having determined the variation and amplitude of the incident 
radiation on the detector cell, the ability of an infrared cell to detect 
this aperiodic signal will be investigated next. A study of response 
characteristics of infrared cells indicated that lead sulphide cells 
have a superior relative response (2 1) in the spectrum of interest 
when compared with cells made up of other photoconductive materials 
such as lead telluride, lead selenide, and indium antimonide. 

Before considering the minimum detectable sional oi a ,celiaam 
will be helpful to define certain terms associated with infrared detector 
application. Noise-equivalent-power, FP, is defined as that incident 
radiation on the cell which produces a signal to noise power ratio of 
unity [7] Another basic parameter, the "Jones quantity," S, is defined 


at some arbitrary reference level by the following relationship, (22) 


= le 
ome)? 


Aln+ 
] 


3 and f,, are the upper and 


lower amplifier bandwidth half-power frequencies, respectively. The 


in which A = sensitive cell area and f 


time constant,t, is defined as the time required for the response to 
a step input to build up to 63 per cent of its maximum value and it has 


been observed that 
St = constant |= 76 


Combining these two equations yields the somewhat simpler and more 


useful expression for noise-equivalent-power, 
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Noise~equivalent~power is generally presented graphically asafunction 
of wave length, as in Fig. 3, for specified reference values of tT, A, 


fand f=. By combining 


Z ] 
P I 2 
. se See (in #) 
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O ] 
O 
and 
ry 2 
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Po = aay) In =- 
t 
where the subscript, "o", applies to the standard reference condition 


and the superscript, "*", identifies those values for a particular cell- 
amplifier application. The noise-equivalent-power for a particular 


xk sk 


*s 3 
cell can be evaluated, when + , A , f, and f, have been deter- 


mined for proper signal output, by the expression 


P + se /2 0 ey 
cor Or O A age | 
Nae ge |) (In £./f,) 
T O pee. ola 
Thus, the noise-equivalent-power, or minimum detectable signal 
power, is directly proportional to the square roots of the cell sensitive 
area and the logarithmic frequency ratio. The quantity is inversely 
proportional to the time constant. In order to evaluate the noise- 
equivalent-power, it will be necessary to investigate the bandwidth 
and time constant required to develop a meaningful signal across the 


detector cell. The majority of military and commercial applications 
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involve modulating the incident radiation as it strikes the detector and 
observing the entire spectrum of the emission source which is visible 
tothe cell. Literature and manufacturers' data are presented with 
these applications in mind and, as a consequence, a special analysis 
will be required to determine the desired cell characteristics for this 
pulse input application. 

The bandpass required to amplify the signal without excessive 
distortion will be determined first. The cell input signal as a function 


of time corresponding to 


Gq == 2000 Btu/ft~ -sec 


has been selected for design purposes. This signal is approximated 


by the parabola shown in Fig. 24 and can be expressed as 


f(t) = 0 20 
f(t) ease > e yee Ai oO” sec 
(a0 eA Geenon meee 


Then, by means of the familiar Fourier integral relations, (23] 


OO 
f(t) = g (w) dw 
- CO 
OO 
gio) = =| £(t)dte 
al © @) 


it is possible to determine the magnitude of the frequency response as 


a function of frequency, from which the required amplifier bandwidth 
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is determined. A minimum bandwidth is desired in order to achieve 


minimum noise-equivalent-power. Proceeding then, 


Oo 40x10~° 
es ~jot | 1 - jut 
g(w) = om it yeclese LS E(t) rdtve 
-cO O 
CO 
ely ~ jut 
oe f(t) dte 
40x10 ~° 
40x107° 
pets - Just 
dee f(t) dte 
O 


Orisa = 40) x oe sec, the approximate test duration in the shock tube, 


we obtain 
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Nondimensionalizing by means of 


cle Te (w)| 5 Jo(w) 


WF O 


3 
lim 2: (3) = $28) (=) 
Ww? Oo 


and expressing lc (0) in decibels permits a simpler form in which to 


where 


evaluate the frequency response for 


] Ww = 10° rad/sec 


which is shown in Fig. 25. The magnitude is down ZO db at 


f, = | Sulke 


which defines a conservative upper frequency limit for the amplifier. 
Representing the signal as a single pulse is somewhat misleading. 
However, an exact solution would have entailed rather tedious com- 
putations not commensurate with the desired accuracy. Inasmuch as 
the system is eventually expected to be used in longer duration shock 


tubes, a wider bandwidth, 


50 x 10° cps 


ry, 
i 


will be utilized. Since the term in the noise-equivalent-power 


expression, ln (ye becomes indeterminate as f, approaches 


] 


zero, the approximation, 


ae a 
fy = ie eee 


will be made for design considerations. 
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The next characteristic to be investigated is the effect of the cell 
time constant on the detector output signal. It should be noted that the 
cell time constant and the noise-equivalent-power are extremely sensi- 


tive to cell temperature, exhibiting the following characteristics:({7] 


Time Constant Noise-Equivalent- Power 
Cell Temperature (microseconds) (watts) 
2526 10 to 1000 Po. say 
(oS 2 te 200 10 Po 
-78°C 100 to 10,000 Po WLS Po 
= NO XS ZUG tomo, cGe AO als P 


A minimum time constant can be achieved by heating the cell, but 
this increases the noise-equivalent-power, reducinp the cellisemcr. 


tivity. The input radiation signal is assumed as before to be 


f (t) 625 oe 0 = treo i0°° sec 


f (t) for all other t 


I 
© 


and a transfer function for the cell is assumed to be of the form, 


] 


7S ees 


where s isthe Laplace variable. Then, using Laplace transforms, 


the output signal can be expressed as 


vow) = {4e5) eee} 
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F(s) = 3 
Ss 
O(s) = H(s) F (s) 
_ 2x 625 1 
a aoe + 1/7) 
= -t/r 
o(t) = £(t)— 12507 [t —r(1 — ) 
The output signal is determined for values of +t = 1, 5, and 10 micro- 


seconds and is plotted in Fig. 26. It is indicated by these curves that 

a minimum time constant is desirable, preferably less than five 
microseconds. This precludes the use of detector cooling and suggests 
heating the cell. However, heating the detector will provide an addi- 
tional source of radiation that will be comparable to that emitted by the 
film, so the alternative of operating with an expected time constant of 
ten microseconds at room temperature must be accepted. Certain 
manufacturers have suggested the possibility of reducing this value of 
+ by a selection process and an attempt to procure a minimum time 
constant cell should be made. 


For this analysis, the time constant will be assumed to be 


+r = 10 microseconds 


The sensitive cell area will be taken as that of the light tube, 


p= Oe > Z oie 


Having determined the three basic parameters, it is now possible to 
evaluate the expected noise-equivalent-power based on the reference 


values in Fig. 3. 


34 





sole x 1/2 
A’ (1n ey) Sp) 
A (1n ele 


(in oe 1/2 
3 19712 (180) (0.1152 1/2 tome 
10 /\ 0.01 ( = 


"995 


CG a Oma ate 


This value indicates that an acceptable signal to noise ratio of the cell 


of six 1s reached at time 


t = 6 microseconds 


when the film temperature 
4 


T = 200ckK 


However, a signal which is distinguishable from the background 


radiation 1s not available until r= 6, Vatewhech time 


If- 


li 380°K 


t = 8 microseconds 


Steady state conditions in the shock tube are reached at about this 
same time, so this means that the period of observation should begin 
at 

t = 10 microseconds 

The magnitude of the output signal from the cell will be deter- 
mined next. This quantity 1s commonly called responsivity and is 
expressed in units of output volts per watt of incident radiation. The 
cell is placed in a series circuit as shown in Fig. 27. Closing and 


opening the switch just before and after each test run will permit use 
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of maximum voltage across the detector without excessive heating of 
the cell. A high bias voltage gives the greatest signal. Hence, an 
applied battery voltage of 400 volts will be used, which is a practical 
upper limit. 

A series resistor is used to match the cellresistance for maximum 
power delivery to the amplifier. The cell resistance varies somewhat 
with manufacture and will be assumed to be one megohm, whichisa 
representative value. The signal from the detector is generated by 
the change in cell conductivity caused by the incident pacrema. 

The responsivity varies with cell area, cell resistance, cell 
voltage, cell temperature, time constant, modulating frequency, and 
source temperature. The responsivity of an individual cell is usually 
determined by experimental means at some standard conditions. The 
data from such tests for a typical lead sulphide cell are presented in 
Fig. 28. The tests consist of focusing the radiation from a constant 
temperature black body source on the cell sensitive area. The tempera- 
ture can be setat any desired level and the black body emitted power is 
distributed over the spectrum in accordance with Planck's spectral 
distribution law. Since the cell only responds to radiation between 
about 1.9 and 3.5 microns, avery small part of the incident energy 
will be in the cell's "visible" spectrum for the lower temperatures. 
However, as the black body temperature is increased, a greater 
portion of the incident power will lie in the detectable wave length 
range, resulting ina greater response for the same total power input. 
For the proposed application in which an interference filter is used, 
all the energy in the selected operating bandwidth lies within the region 


of maximum response of the detector and hence all filtered incident 
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radiation will be "seen" by the cell. It follows, then, that the cell 
responsivity should correspond to the maximum temperature line on 
Fig. 28. The maximum response for the cell closely approximates 
the 2000-3000°K black body curve, and this will be used for design. 
The responsivity given by this curve must be corrected for those 
parameters differing from the standard conditions, for example, the 
time constant, cell area, and cell resistance. Typical correction 
curves are shown in Figs. 29 and 30. Evaluating the following 


expression, 


Ridb) = Retandard OS A sae const f ecou area resistance 


yields a cell responsivity of 
R(db) = 404+ 15+5 = 60 db/microwatt 
which is 
R = 1000 volts/watt 


For the selected design radiation signal of Fig. 24, a signal level of 
approximately 100 to 1000 microvolts is predicted for the detector 
output. 

This signal must be amplified and transmitted to the end device. 
The amplifier should have a bandwidth from 0-50 kilocycles and 
amplify the microvolt signal to an acceptable level for coaxial cable 
transmission. An oscilloscope, such as the Tektronix 535, should 
serve quite satisfactorily as an end device, with the addition of a 
Land Camera, to record the output trace. 

A most critical problem will be the prevention of stray pickup 


by the transmission leads and by the circuitry within the shock tube 
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instrument housing. Radiation from the stray electromagnetic fields 
will be quite extensive. The 25-kilovolt discharge in the shock tube 
and the intense magnetic field from the nose coil and the associated 
leads will generate a great deal of this spurious radiation. To shield 
the oscilloscope, it willbe placed in a copper screen room about 

20 feet from the shock tube. The detector-amplifier unit will be 
enclosed in the copper housing which leaves only the coaxial cable 
lead from the housing as a source of stray pickup. The use ofa 
differential amplifier at the oscilloscope will minimize the effect 

of this stray pickup. Elementary tests should be conducted to 
determine whether or not additional shielding is required. 

The amplifier unit should have a high input impedance matching 
the detector cell resistance. An amplifier gain of approximately 
1500 should be ample. Power supply oe is extremely critical 
and the undesirability of operating the amplifier with a remote 
supply, due to stray pickup, suggests the use of transistors. The 
direct-coupled amplifier design in Fig. 31 with a grounded collector 


input stage should prove to be adequate. 
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7. Calibration. 

The analysis of the preceding sections suggests that the magnitude 
accuracy of the measured heat transfer rate will be somewhat in doubt. 
The transmission path attenuation values, expectedcellandinterference 
filter performance and many other factors were estimated as closely as 
possible, but to measure the heat transfer rate quantitatively, it 
becomes necessary to devise a means of calibrating the system. The 
most obvious method is to subject the film to an electrical energy input 
that dynamically simulates the heat energy input of the shock tube tests, 
and record the detection system response. The joule heating of the 
gold film will simulate the heat flux and it should be possible to record 
various readings for various magnitudes and durations of current pulses 
thaouch the film. [The peak current required jor me oueu A thick gold 
film appears excessive for usual methods of current pulsing, due to the 
low resistance offered by the film. Extrapolating the findings of 
Holland, [2.4] it is expected that for film thicknesses somewhat greater 
than 1000 A, the bulk resistivity for gold is valid. From the resistance 


relationship, 


R = p 4 ohms 


and assuming a square film, 0.38 centimeter ona side, the resistance 
Po 1OUnd tO be approxinately 5 x 10° ohms. 

Simulating a heat transfer rate of 5000 Btu/sec-ft- would require 
650 watts of power or about 115 amperes discharged through the film 
for the desired time interval. The ditticulty in controlling this mmacna. 
tude of current gives rise to the undesirable, but perhaps necessary, 
alternative of calibrating the system with a thinner film, and then 


replacing the 5000 A film for measurements. High currents of long 
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time duration would very likely physically destroy the film, or arc 
over at the contacts, preventing dissipation of the desired energy input. 
Hot spots would probably disturb the homogeneity of the film surface 
and destroy the infrared opacity. Based on values reported by Holland, 
a 100 A thick gold film will have a resistance of approximately 40 ohms. 
Calibrating with this film will permit dynamic heating of a foil to the 
expected experimental temperatures with a current of about three to 
four amperes. The radiation given off by the 100 A film should closely 
approximate that energy emitted by the 5000 A film. Interchanging the 
film-substrate configuration for calibration purposes will undoubtedly 
lead to a variation in transmission path efficiency, due to physical 
differences in the interfaces between the gold film and the substrate, 
and between the substrate and the sapphire light tube. Also the physi- 
cal properties of the quartz substrates may not be identical. 

However, a calibration would seem of interest even if the accuracy 
is only within 50 per cent of the true value. The 100 A film can be 
pulsed for a desired time interval with various current magnitudes 
corresponding to different heat transfer rates. Equating the joule 


heating rate to the heat transfer rate, 


the magnitude of the current pulse in terms of time and heat transfer 


rate can be determined by 


iicpet lta amperes 
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where 
= Ro (ea aes) 
a = thermal coefficient of 
resistivity 


It is desired to express the change in temperature in the expression 
above as a function of time. For very thin films, as shown in Sect. 3 
and elsewhere, [25] the temperature change of the film can be expressed 


as, 


Q(x, t) = AT (x, t) =a Tt for | aeO 
rc, p 


where the physical properties are those of quartz. 


AT(0,t) = 8.5 qt (°K) 


\x | : 
Me) Cots) a = (1 + 8.5 aq|t) 1/2 amperes 
O 


Assuming the bulk value for the thermal coefficient of resistivity is 


then 


valid, the current pulses corresponding to different heat transfer rates 
for a film resistance of 100 ohms are plotted in Fig. 32. Thisindicates 
that a current supply network capable of one to five ampere pulses, 
with a duration of 15 to 40 microseconds, should prove adequate for 
heat transfer rates of interest. Average current values will be used 
here approximating the actual curves of bie. 

In order to handle the magnitude of currents required in the 
desired time duration, the use of electronic switching is suggested. 
Thyratrons with the required current carrying capacity such as 


type 3D2Z2A can be utilized as switches. In order to supply this 
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current, conventional power supplies of milllampere capacity will not 
suffice. It seems logical to store the required charge on a high voltage 
capacitor. The circuit shown in Fig. 33 was suggested by the authors 
and appears to meet the requirements. The two thyrairens are mor. 
mally off, with the control grids held negative by the grid bias. The 
Capacitor, Co is charged to the desired plate voltage, with the 
switch in position No. 1. Whenthe capacitor has been charged, the 
switch is shifted to position No. 2, the desired current magnitude set 
up by means of the potentiometer, Bae and the pulse duration, At, 
selected on the pulse generator and delay network. When the pulse 
eenerator 18 triggered, pulse No. | turns thyratron No. ion pen. 
Imitting current to flow through the gold film. Pulse No we2 tum. 
thyratron No. 2 on after the desired pulse duration, and shorts the 
PeLrenttometer and cold mlmuronc momma. 

However, experimental determination of the conductivity charac- 
teristics of the 100 A thick film should be made prior to actually 
pulsing the film. The conductivity of the film should be investigated 
over the temperature range of interest, and the values obtained used 
to determine current pulse characteristics such as those shown in 
Fig. 32. When this has been accomplished, then the actual calibration 
can take place, utilizing the described network. 

Other means of calibration that could be investigated are: 

(1) Heating a film on a separate substrate, and then placing the 

configuration in contact with the sapphire rod 

(2) Alternating current joule heating 


(3) Ignitron current pulsing. 
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Suggestions 1 and 2 are static methods, but would seem simple to 
perform, and could give quite accurate calibration for known film 
temperatures. A thermocouple could be utilized to provide precise 
film temperature measurements at the time of calibration. 

Method 3 offers the hope of dynamic calibration of the 5000 A 
film with electronic control. However, as stated previously, whether 
or not the film can physically withstand current pulses of the required 


magnitude remains to be determined. 
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8. Conclusions. 

It is concluded from this analysis that the proposed infrared 
system will be capable of measuring stagnation point heat transfer 
rates on test bodies subjected to simulated space flight conditions 
in an electromagnetic shock tube. The ultimate design consists of 
a thin metallic film evaporated on a quartz window. Radiation 
through the quartz window is transmitted via a sapphire light pipe 
to a sensitive lead sulphide detector cell. The output of the cellis 
amplified and observed on a conventional oscilloscope. 

It was found that optimum system performance will be achieved by 

(1) Use of a gold film evaporated on a fused quartz substrate 

(2) A film thickness of 5000 A 

(3) A sapphire light pipe as a most efficient means of energy 

transmission in the near infrared spectrum 

(4) Selection of a lead sulphide cell 

(5) Choosing a spectral bandwidth of 0.4 micron centered ata 

wave length of approximately 2.8 microns 

(6) An interference filter to limit the lower wave length energy 

striking the detector cell 

(7) Electromagnetic shielding of all electrical circuitry. 

The apparatus will be useable under the ionized flow conditions 
in which present resistance film gages are inoperative. It is believed 
that heat transfer rates of approximately 2000 to 35,000 Btu/ft“-sec 
can be measured within the 15 to 40 microsecond observation periods. 
The upper value is limited by the melting temperature of the gold film 


or the detector cell time constant. The lower limit is the minimum 
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heat transfer rate which will produce a detectable film temperature 
within the maximum available test time. 

The accuracy of the absolute heat transfer rate measurement will 
be somewhat doubtful, until a more appropriate means of calibration is 
devised. The proposed method should suffice for preliminary experi- 
mentation and relative measurements, but more precise calibration of 
the apparatus will be of future value. Further calibration studies 
should be undertaken to achieve the ultimate system capability. 

Investigation of the following subjects of current interest could be 
undertaken with this apparatus. 

(1) Radiative, convective, and recombination heat transfer 

to bodies in hypersonic flow 

(2) Effect of magnetic fields on radiative, convective, and 

recombination heat transfer rates. 

The system will be assembled and ready for experimental work 
during the summer of 1959. Experimental difficulties remain to be 
found and remedied, but it is believed no major problem should arise. 
The first application of the apparatus will be to verify the theoretical 
work done in magnetoaerodynamic cooling. The successful completion 
of this experimental work should greatly contribute to future space 


vehicle design and performance. 
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Figure 2. Details of the Shock Tube with the Test 
Model Installed. 
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SHOCK WAVE 





Stationary Model Coordinate System 





Moving Model Coordinate System 


Figure 5. Shock Wave Coordinate Relationships 
within the Shock Tube. 


a4 





SHOCK WAVE 





Stationary Model Coordinate System 





Moving Model Coordinate System 


Figure 5. Shock Wave Coordinate Relationships 
within the Shock Tube. 


ye 





P, /Po 


FREE STREAM DENSITY RATIO, 


002 


O 


Figure 6. 


008 ,40000 
006 A 30000 


004 





SHOCK MACH NUMBER 


Free Stream Density Ratio and Free Stream 
Velocity in the Shock Tube. 


310, 


FT/SEC. 


FREE STREAM VELOCITY, 





HEAT TRANSFER RATE, BTU/FT? -SEC 


16000 
12000 Ee 
8000 BY 


10. 20 30 40 50 
SHOCK MACH NUMBER 


Figure 7. Predicted Heat Transfer Rates. 
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Figure 9. Ratio of Emitted Energy to Transmitted Energy, 
W/ We for Various Film Thicknesses. 
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Figure 10. Heat Flow in a Semi-Infinite Slab. 
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Temperature-Time Profile in a Semi-Infinite Slab. 
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Figure 12. Outer Surface Temperature for a Semi-Infinite 
Slab as a Function of Time. 
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Figure 13. Heat Flow Through a Parallel Wall Slab. 
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Figure 14. Interface Temperature-Time Relationships. 
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Figure 16. Sapphire Spectral Absorption 
Coefficient. 
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Figure 17. Transmission Efficiency of a 
Sapphire Light Pipe. 
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Figure 18. Detailed System Assembly. 
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Figure 19. Eddy Current Heating Considerations. 
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Ratio of Emitted Power to Extraneous Power 
Incident on the Detector Cell. 
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Figure 21. Desired Interference Filter Characteristics. 
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Figure 22, Ratio of Emitted Power 1o Pxtraneous over 
Incident on the Detector Cell at a Wave Length 
of 2.56 Microns. 
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Figure 23. Incident Signal Power on the Detector Cell 
as a Function of Time and Heat Transfer 
ater 
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Figure 24. Approximating the Incident Power Relationship. 


70 





OOOOO! 


‘sjusursrimbey easuodsey Aouonbea gq szatsttduy 


OOOO! O00! OO! Ol 
‘95S Y3d SSIDAD ‘AON3NO|3YS 


"GZ earns gy 





“TWNOIS 


S17381930 


(pe 





POWER, MICROWATTS 


O 


a0 


oO 


“CT 
Bees, 
aan 
me 


fs 


eZee 





ATL 
—_ 


O 1O 
TIME, MICROSECONDS 


Figure 26. Detector Cell Time Constant Response. 
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Figure 27. Two Possible Detector Cell Circuits. 
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Figure 28. Typical Lead Sulphide Cell Responsivity Characteristics. 
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Figure 29. Correction to Responsivity Signal for Modulating 
Frequency and Time Constant. 
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Figure 30. Correction to Responsivity Signal for Cell Area 
and Cell Resistance. 
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Figure 3l. 


Direct-Coupled Transistor Amplifier. 


77 





AMPERES 


CURRENT, 


4.0 


3.0 


—_ Tr 


PAF An 





TIME, MICROSECONDS 


Figure 32. Calibration Current Pulse Requirements. 


78 













PULSE 
GENERATOR 
AND 
TIME DELAY 


aN’ 3D022A 
Ra, 
ne H.V. 
ATH } C, POWER 
SUPPLY 
Ra, , 
= 


3D22A 
PULSE 7| 7 


o 


ae Aa 


PULSE "2 | 


’ 
@] { —— 









PULSE ”2 


















RADIANT 
eee 
ENERGY 







IR-DETECTOR|,$ 


7 


OSCILLOSCOPE 


Figure 33. Proposed Calibration Circuit. 
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